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ABSTRACT 

Cosmic rays and magnetic fields can substantially impact the launching of large-scale galactic winds. Many 
researchers have investigated the role of cosmic rays; our group previously showed that a cosmic -ray and 
thermally-driven wind could explain soft X-ray emission towards the center of the Galaxy. In this paper, we cal- 
culate the synchrotron emission from our original wind model and compare it to observations; the synchrotron 
data shows that earlier assumptions about the launching conditions of the wind must be changed: we are re- 
quired to improve that earlier model by restricting the launching region to the domain of the inner "Molecular 
Ring", and by decreasing the magnetic field strength from the previously assumed maximum strength. With 
these physically-motived modifications, we find that a wind model can fit both the radio synchrotron and the 
X-ray emission, although that model is required to have a higher gas pressure and density than the previous 
model in order to reproduce the observed X-ray emission measure within the smaller 'footprint'. The drop in 
magnetic field also decreases the effect of cosmic -ray heating, requiring a higher temperature at the base of the 
wind than the previous model. 

Subject headings: ISM:outflows - ISMxosmic rays - ISM:magnetic fields - Galaxy:kinematics - 
Galaxy:evolution - Xrays:diffuse background 



1. INTRODUCTION 

The possibility of thermal-pressure driven outflows on 
galactic scales has b een considered for some time (see, e.g., 
IVeilleux et al.ll20 05, and references therein). Somewhat less 
well-known, however, is the possibility that cosmic rays and 
magnetic fields can also help drive outflows: if, as commonly 
assumed and even sometimes observed, magnetic-field pres- 
sure and cosmic -ray pressure are in approximate equipartition 
with the hot-gas thermal pressure within galactic disks (e.g., 



Duridl990tlP^[T^lHeileslp996l lZ weibel & Heilesl ll997t 
Webb erl 1 998t lBeckll200 It lCoxfl2005l) . these sources of energy 
will be important in large-scale outflows in more quiescent 
galaxies, such as perhaps our Milky Way. Understanding and 
constraining the dynamics and observational constraints on 
such kpc-scale galactic outflows is the main aim of the present 
paper. 

The possibility of a kpc-scale outflow from the Milky 
Way h a s bee n studied for over three decade s (e.g. , 
Ipavi chl 119751; IBreitschwerdt et al.1 119871 11991 " 
Bland-Hawthorn & Cohenl 120031: iBreitschwerdtl 12003 



Veilleu x et al. 1 120051) . including work on the observational 
effect s of such an out fl ow (e .g. . iLerche & Schlickeiserl 
1981 Uokipii & MorfiUl [19871 IReich & Reich! ll988aT 
Pohl & Schlickeiserl Il990t iBloemenl 119911; iBloemen et al l 
1993; Breitschwerdt & Schmutzler 1994; Zirakashvili et al. 
1996t iPtuskin et al.1 119971; IBreitschwerdt et al"n 2002). 
Bland-Hawthorn & Cohenl (120031) showed possible evidence 
for such an outflow in X-ray, infrared, and radio observations, 
along with models of a bipolar structure. In Everett et al.1 
(2008, hereafter Paper I), we continued this work by building 
a wind model powered by both cosmic -ray pressure and 
thermal-gas pressure. In such winds, cosmic rays yield 
some of their energy to Alfven waves on the magnetic field 
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via the "Streaming Instability" ( Kul srud & Pearcel 119691) 
which are then damped to deliver energy and momentum 
to the thermal gas. We applied that model to observations 
of excess diffuse X-ray emission towards the in ner Galaxy 
dSnowden et all [19951 \l99H. lAlmv et alj p000). Paper 
built especially on IBreitschwerdt et al J ( 1991 ). 1 
the following modifications: the wind is launched from a 
restricted range of Galactocentric radii (an annulus in the 
inner Galaxy), originates at the Galactic midplane with a 
vertical magnetic field, and (as indicated by our analytical 
work) the cosmic -ray generated Alfven waves are completely 
damped, subsequently heating the gas. This model fit the 
longitude-averaged X-ray emission observed with ROSAT 
near 0.65keV and 0.85keV with \ 2 val ues a factor of tw o 
smaller than the previous static model (lAlmy et al.l [2000), 
while simultaneously yielding pressure and density parame- 
ters approximately equal to those already infe rred for hot gas 
near this region in the inner Milky Way (c .f . . iFerrierd 1200 ll 
hereafter, F01). 

The above initial model (reviewed in © was derived solely 
from a fit to the soft X-ray emission and made several as- 
sumptions that bear continued examination. In this paper, 
we test the wind model further by comparing calculations of 
the synchrotron emission from the wind to 408-MHz all-sky 
surveys towards the center of the Galaxy: radio-synchrotron 
observations are an important constraint on a cosmic -ray 
and thermal-gas pressure-driven outflow, given the depen- 
dence on both the cosmic-ray pressure (albeit through the un- 
certain and energy-dependent ratio of cosmic -ray electrons 
to protons) and the magnetic field strength, both of which 
are crucial in our model. Many groups have worked to 
understand the distribution of synchrotron emission in t he 
halo (TO 1 : iBeuermann et al.l 1 1 9851: IReich & ReicE 1 1 988allbl : 
Broadbent et alj|1990t IStrong et aUboOOt IReich et al.ll20o4 
Sun et al.ll2008Ude Oliveira-Costa et al Jl2008t IReich & Reich! 
20081; [Waelkens et al. 2009), but more work is needed to un- 



derstand the role of cosmic-ray advection into the halo via a 
wind and the subsequent cooling of cosmic-ray electrons via 
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inverse-Compton and synchrotron cooling (although we note 
that the cosmic-ray propagation code GALPROP has an ap- 
proximate Galaxy -wide wind model with a constant velocity 
gradient; see, e.g.. lPorter et alJl2008l and references therein). 

In this work, we build a model of the synchrotron emis- 
sion of a large-scale wind launched from an annulus in the 
inner Galaxy, taking into account synchrotron and inverse- 
Compton cooling and the full velocity gradient in the wind; 
this work is detailed in $3] We also present predictions for the 
synchrotron spectral index, but note the possibly strong im- 
portance of foreground synchrotron emission. By comparing 
the observed synchrotron emission to that calculated for the 
model of Paper I, the observations can then be used to exam- 
ine two assumptions made in Paper I: the radial range on the 
Galactic disk where the wind is launched and the magnetic 
field strength. 

The calculated synchrotron emission will show that the 
wind model of Paper I overpredicts the observed radio emis- 
sion; the observed synchrotron emission is therefore a strong 
constraint on the wind model. As outlined above, we then 
consider if it is possible for a modified wind model to fit the 
observations; these possibilities are explained in §|4] First, the 
wind may have a thin ner footprint than the AR = 3 kpc width 
used in Paper I ( 34. Il l; this is suggested b y observations of the 
"5 kpc Molecular Ring" (Jackson et al. 2006, B. Benjamin, 
private communication). Second, the magnetic field strength 
within the wind might be lower than the value of 7.8//G used 
in Paper I ( 34.21 i. Third, the cosmic-ray electron-to-proton ra- 
tio may be different towards the center of the Galaxy, although 
we find it unnecessary to invoke this, given other, more physi- 
cally motivated and astrophysically constrained modifications 
to the model ( 33. U . (We also modify the x 2 fit parameter to 
avoid biases du e to assumptions about X-ray absorption in 
34.31 ) Then, in 34.41 we present the best-fit wind model that 
addresses both the soft X-ray and synchrotron observations, 
and compare this model to that of Paper I. Finally, we discuss 
the results in $5] 

2. THE HYBRID COSMIC-RAY AND 
THERMALLY-DRIVEN WIND MODEL 

Before we outline our model of synchrotron emission from 
the wind, we first review the wind model from Paper I. This 
model was motivated by observations of the inner Milky Way, 
where an excess of T ~ 3 x 10 6 K, X-ray e mitting gas has 
been observed with a scaleheight of ~2 kpc (ISnowden et all 
fl99l fl997t iPark et al.l[l997l fl99l lAlmv et al.ll2000h. Some 
attempts had been made (Snowden et al. 1997; Almv et al. 
2000) to fit this emission with static gas distributions, as 
the gas pressure of the plasma is insufficient, on its own, 
to drive a wind from the Milky Way. Our group asked 
whether cosmic-ray momentum and energy could be com- 
municated, via the "Streaming Instability" (IWentzell 119681 : 
iKulsrud & Pearcdll969t IKulsrud & Cesarskyi U97flT to the 
thermal gas, and hence be used to launch a wind. We 
were not the first to ask this question; building on much 
past research into the power of cosmic rays to affect 
gas dynamics and help drive a wind (e.g.. Ilpavichl [19751: 
Skilli ngl 11975b iBloemenl 11991b iB reitschw erdt et al.1 ll99lL 
1993; Breitschwerdt & Schmutzler 1994; Zirakashvili et al. 
19961 iPtuskin et al.lll997HZirakashvili & Volkll2006l) . we as- 
sembled a ID hydrodynamical model of a cosmic -ray pres- 
sure and thermal-pressure driven wind. The equations for this 
wind model are reviewed in the Appendix. 



The model of Paper I differed from previous models in 
that we concentrated solely on understanding the enhanced 
X-ray emission towards the inner Milky Way. Taking the hy- 
pot hesis that the em ission is centered on the Galactic Cen- 
ter (lAlmy etalJl2"OO0h . we limited the maximum radial ex- 
tent of the wind to 4.5 kpc in Galactocentric radius. We 
assumed the wind was launched from the midplane of the 
Galaxy; the wind then flows vertically within flowtubes with 
cross-sectional area given by a simple analytic form with a 
characteristic expansion h eight that is a free parameter (as in 
IBreitschwerdt et aP 19911) . We set the initial cosmic -ray pres- 
sure and magnetic field strength to the values inferred from 
synchrotron observations at the Galaxy's midplane (F01), as- 
suming a vertical magnetic field. The initial gas pressure 
and density at the Galactic midplane, which are relatively 
poorly known, were left as free parameters. In early tests, we 
found that for fiducial values of all these parameters, winds 
could not escape from the very inner portion of the Galaxy 
(R < 1.5 kpc), so we hypothesized an annular geometry, with 
the wind only occupying a thick ring from 1.5 kpc to 4.5 kpc. 
(This inner region without a wind may correspond to the tran- 
sition from a wind to an X-ray emitting atmosphere of gravita- 
tionally confined sup ernova-heated gas hypothesized in more 
massive spheroids by lBinneyll2009l) . 

Values for the three free parameters (gas pressure & density 
at the midplane, and the scale-height of flowtube expansion) 
were then found by fitting the predicted X-ray emission of the 
wind to the observed emission in ROSAT's R4 and R5 bands; 
the wind was found to provide a significantly impr oved fit (by 
a facto r of two in \ 2 ) over the static poly trope of Almv et al. 
(2000). In addition, the gas pressure and density values that 
were found in this process were plausible for the inner Galaxy 
(see TableQ]). Finally, the power required to launch this wind 
was approximately a factor of two higher than the estimated 
supernova power, which was plausible given the many ap- 
proximations of the theory, but possibly indicated that further 
model refinement and/or constraints were necessary. Despite 
the clear simplicity of this approximate model, such a cosmic- 
ray pressure- and thermal-gas pressure-driven wind within our 
Galaxy seemed to fit the observations well. 

X-ray observations alone do not uniquely determine the 
structure of the wind; comparisons of the model to syn- 
chrotron data allow a natural and complimentary check. We 
outline the synchrotron-emission calculation and compare the 
predicted emission to the observations in $3] The predicted 
synchrotron emission is larger than what is observed, so we 
vary other aspects of the wind model in $4] The resulting 
best-fit model provides nearly as good a fit to the X-ray obser- 
vations as the original model of Paper I and is more consistent 
with other known aspects of galactic structure. 

3. CALCULATING THE SYNCHROTRON EMISSION 

Given the above magnetic field and cosmic-ray density (de- 
fined by the wind model presented in Paper I), the synchrotron 
intensity al ong a line of sight can be integ rated via equation 
(3.20) from lGinzburg & SyrovatskiH (119651) : 

V3e 3 f ( ' v f°° \ 
I v = T / n(£e,r,k)B(r)sin6» — / K 5 Mr))dr) dE e dr. 

m eC Z J \VeJvlv, J 

(1) 

where e is the electron charge, m e is the electron mass, and 
c is the speed of light. n(E e ,r,k) is the cosmic -ray electron 
density, which is a function of cosmic-ray electron energy, 
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E e , position, r, and direction, k. ^5/3 is the Bessel function 
of the second kind, B(r) is the magnetic field strength as a 
function of position, 9 is the angle between B and k, v is the 
observation frequency, and v c is defined by: 



3eB_ 
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vert this result from the Sun's position to a Galactocentric 
radius of R = 3.5kpc by multiplying that local cosmic-ray 
electron flux by a factor of 2.4, as derived in Equation 11 
of F01. This factor comes from fitting the variation of syn- 
chrotron radiation in th e Galactic midplane as observed by 
iBeuermann et al.l (11985b . 

For consistency, we also use the cosmic-ray electron- 
spectrum power-law of 7 = 2.5 from F01. We assume that 
this cosmic-ray electron spectrum is constant at the Galactic 
midplane, and constant over time, evolving only with height 
in the wind due to cooling. Of course, if the CR e"/p ratio 
changes towards the center of the Galaxy from that observed 
locally, such a change would produce a directly proportional 
change in the synchrotron intensity. We point out that a differ- 
ent input spectrum may be possible, as IStrong et al.1 |2004a) 
find that their Galactic propagation models imply a different 
electron-injection spectrum from that used in Paper I or in the 
present work. Such uncertainties represent a large area of pa- 
rameter space that adds more complexity via another free pa- 
rameter, and so we choose instead to retain the present spec- 
trum, based on local measurements. Our synchrotron results 
could easily be scaled to lower or higher CR e" densities; as 
we will see in §|4] this is not necessary to fit the observed syn- 
chrotron emission. 

We note briefly that the above procedure is similar to that 
we used to estimate the cosmic -ray proton pressure at the base 
of the wind in Paper I; as in F01, we have assumed that the 
local cosmic -ray proton to cosmic -ray electron ratio applies 
also to the are a at the base of this wind. This is mentioned 
further in $4.41 The electron-to-proton ratio is set to 1 :50; one 
could assume various ratios up to 1:100, perhaps, but assum- 
ing 1 :50 should force relatively strong constraints on the wind 
model when comparing to radio synchrotron emission. 

3.2. Cooling the Cosmic-Ray Electrons 



To include the effects of inverse-Compton cooling 
on the cosmic-ray electrons, we use the representation 
of the interstellar radiation field (ISRF) that has been 
developed for use with the GALPROP code (avail- 
able at http://qalpr o p.stanfor d .edu / ; see 
Strong & M oskalenko 1998; StrongetaD l2000l l2004at 



1 ; IPorter & et alj|2005t iPorter et a l. 20080 In their model, the 



For evaluating and integrating the full synchrotron intensity 
equation, n(J? e ,r,k) is required throughout the wind: we cal- 
culate »(E e,r,k ) by assuming a given initial electron dis- 
tribution Q3.U . and accounting for the int egra ted inverse- 
Compton, synchrotron, and adiabatic losses ( $3.2t as the wind 
is launched from the disk. After we give the details for the 
calculation, below, w e br iefly outline the converge nce tests 
the code has p assed ( 93 .3t . define the backgrounds ( 93 .4t and 
foregrounds ( $3.5t that must be included, and sho w th e pre- 
diction for synchrotron emission from the wind in $3.61 

3.1. The Galactic Midplane Cosmic-Ray Electron Spectrum 

The wind equations for this model consider both the ther- 
mal gas and cosmic rays as fluids. To model the synchrotron 
radiation of the cosmic-ray electrons, we have to model in 
more detail the cosmic-ray electron spectrum with energy. To 
do this, we first take, as a typical cosmic -ray electron spec- 
trum, the electron spectrum observed at the Galactic midplane 
at the Sun's position (F0 1). As i n F01, we use the normal- 
ization from Figure 6 of I Webber! 0983) where, at 2.3 GeV, 
dn/dE = 17.3m" 2 sec" 1 ster" 1 GeV/nucleon~ . We then con- where 



ISRF is given as a function of wavelength on a regularly 
spaced grid in Galactocentric radius and height above the 
plane: the radial grid has 41 elements in Galactocentric 
radius, starting at 0.25 kpc and spaced at 0.5 kpc increments. 
The height grid starts at z = -5 kpc and is spaced in incre- 
ments of 0. 1 kpc, in 101 steps, to z = 5 kpc. The ISRF photon 
field is then converted to f~, {uS), the background photon 
distribution (see Appendix C, [Strong & Moskal enkoTl 1998t) 



via the definition: 



t/ph = m e c' 



00 f 7 (uj)duj , 



(3) 



where t/ pn is the energy density of the ISRF, uj is the pho- 
ton frequency, and 7 is the electron Lorentz factor. With 
the ISRF photon distribution, we can then write the equa- 
tion for inverse-Compton co oling (see, again, Appendix C of 
Strong & Moskalenko 1998, and references therein) as: 



dE 
dt 



ic 
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(4) 
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(5) 



and where £_(•••) denotes t he dilogarithmic function (th e 
polylogarithm with n=2; see Abram owitz & Stegunl 119721) . 
f3 = v/c for the cosmic-ray electrons, and fc^ = ui'yil ± /?). For 
completeness, we note that the full application of this formula 
to the production of 7-ray e mission requires the considera tion 
of anisotropic scattering (see Moskalenko & Strong 2003); in 
the present case of simple cooling, where we need only the 
absolute power emitted, the anisotropic nature of the radia- 
tion field is not important. 
We also include the energy loss from synchrotron cooling: 



dE\ 
It) 



32 



7rr e z c£/ B 7 2 /3 2 - 



(6) 



synch 



Adiabatic cooling is already included via the hydrodynamic 
cosmic -ray model of the wind. 

In summary, then, the hydrodynamic wind model in- 
cludes the adiabatic cooling, and then a separate synchrotron- 
modeling program implements the energy-dependent inverse- 

4 The ISRF data file is http : //galprop . Stanford. edu/FITS/ 
MilkyWay_DRO . 5_DZ0 . 1_DPHI10_RMAX2 0_ZMAX5 
_galprop_f ormat . f its . gz, downloaded on September 23rd, 
2008 from http://galprop.stanford.edu/web_galprop/ 
galprop_home . html. 
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Compton (Eq. |4|i and synchrotron (Eq. |6]l cooling losses for 
the energy resolved electron distribution. 

To calculate the effect of these cooling terms, the At for 
each discrete vertical step of the wind is calculated from the 
velocity profile, and then simply multiplied by the above 
cooling rate at every step. The inverse-Compton and syn- 
chrotron cooling do substantially impact the energy density 
of the cosmic-ray electrons. It is important to note, however, 
that neither process impacts the cosmic -ray protons, which 
are the source of momentum for driving winds; cosmic-ray 
electrons, with a number density only approximately 2% that 
of the cosmic-ray protons, do not contribute significantly to 
the cosmic -ray pressure. So, in this case, when we consider 
cosmic-ray electrons as a separate component, and cool those 
electrons in an already defined wind model, that cooling does 
not impact the wind model (such as the mass outflow rate and 
velocities achieved, for instance). 

Following the calculation of the local electron spectrum, 
the local intensity is calculated and then integrated numer- 
ically on lines of sight through the wind; synchrotron self- 
absorption does not play a role in the Galaxy for v > 50 MHz, 
so it is neglected here. This calculation is carried out over 
a grid in Galactic longitude & latitude. We then average the 
emission over longitude (as done previously for the compari- 
son with the soft X-ray emission) and plot it as a function of 
latitude to compare with observations. 

How do the various cooling terms compare? We have found 
that inverse-Compton and synchrotron cooling of cosmic-ray 
electrons dominate adiabatic cooling by a factor or ~ 10 
near z ~ 500pc, but that adiabatic cooling achieves parity 
around z ~ 2kpc above the Galactic midplane, reaching a 
maximum in relative strength (2.5 times the total of inverse- 
Compton and synchrotron cooling) at approximately 5kpc; 
above this height, adiabatic cooling decreases again to be- 
come approximately equal to both inverse-Compton and syn- 
chrotron cooling (we have calculated the above energy losses 
over the entire electron-energy distribution that we consider, 
from ~ 10~ 3 GeV to 100 GeV). The models also demon- 
strate that inverse-Compton cooling is more important than 
synchrotron cooling by approximately an order of magnitude 
near R = 3.5 kpc. Also, we have found that Klein-Nishina ef- 
fects are significant at 408 MHz, with total inverse-Compton 
cooling dropping to 65% of that calculated from the non- 
relativistic approximation. 

3.3. Testing the Synchrotron Code 

Before running the full synchrotron calculation with cool- 
ing, we developed an initial version of the code that assumed 
a constant (non-cooling) cosmic -ray electron population with 
a power-law index of 2.7. This code was checked against an- 
alytic estimates for the synchrotron emission along various 
lines of sight, and was used to check convergence of the flux 
integration with various line-of-sight stepsizes; the algorithm 
reproduced the analytical results within ~ 5%. 

We have also tested the synchrotron code by checking for 
convergence in a wide range of parameters used within the 
code: for instance, we have checked that the code converges 
when using a larger number of electron bins in the CR- 
electron spectrum, smaller vertical step-sizes in the cooling 
calculation, and higher resolution in all numerical integrals. 
All of the higher-resolution studies yield less than 5% differ- 
ences in the predicted synchrotron emiss ion. This algori thm 
also reproduces (within 3%) the result of Ferriere (1200 lb for 
the synchrotron emission of the Galactic midplane. Finally, 



the code can also be run in a "static halo" mode without any 
cooling terms or flow velocities; those runs reproduced the 
expected flat spectral indices to approximately 1 part in 10 4 . 

3.4. Synchrotron Backgrounds 

Before comparing to observations, we must also account 
for the various backgrounds in the data, adding those compo- 
ne nts to the w i nd mo del's synchrotron emiss ion. As shown 
in IReich et al.l (120041) . and in more detail in iReich & Reichl 
(2008), the observed brightness temperature (7obs) can be un- 
derstand as the sum: 

Toes = ?gal + Tcmb + Texg + ?off (7) 

where 

Tqal = intrinsic Galactic brightness temperature (8) 

?cmb = cosmic microwave background = 2.73 K (9) 

Toff = zero - level error (10) 

Texg = unresolved extragalactic sources (11) 

We calculate the extr agalactic contribution from the formula 
(see IReich et alJl2()04l and references therein): 

^ =30K (l 78MHz)' 2 ' 9 (12) 
and take Toff from lReich et alj (|2004). 

3.5. A Strong Local Synchrotron Foreground? 

After adding in all of these background components to the 
model, we were intrigued to find a relatively constant level 
of high-latitude (at \b\ > 40°) emission that remained un- 
explained, and would not be explained by the wind emis- 
sion, which would only be im portant up to \b\ y 20° . The 
Galactic synchrotron model of Beuerm ann et alJ (119851) also 
does not address this component; they only fit the syn- 
chrotron emission at low latitudes (b < ±3°). So, it seems 
that the nature and origin of this high-latitude component 
are not well understood, a nd have not been accou nted for 
in the "standard model" of Beuermann et aTl (119851) . From 
the near constancy of the high-latitude emission, from the 
recent evide nce for a relatively local synchrotron emission 
component jFle ishman & Tokarevl 119951; iRoger etafl Il999t 
IWollebenl 120071; ISun et alj 120081 and W. Reich (personal 
communication)), and from our own (brief) analysis of the 
408 MHz to 1420 MHz spectral indices which seem to indi- 
cate a temperature spectral index near j3 ~ 2.75, we conserva- 
tively hypothesize that this high-latitude emission at \b\ > 40° 
is local, and therefore is a foreground Galactic component 
that should be accounted for separately when comparing the 
wind to the observed emission. This is important, as the wind 
model should not predict a larger synchrotron emission than 
is observed in any part of the sky; the presence of a local syn- 
chrotron foreground makes this, of course, more difficult for 
the model, which is therefore even more restricted by the ob- 
servations. 

In addition, as we will show later, these foreground compo- 
nents are also essential in comparing the synchrotron spectral 
index of models to observations, and so highlight the impor- 
tance for a more global model of the Galactic synchrotron 
emission (perhaps like that i n GALPROR but including a 
more general wind model; s ee IStrong et"aTll2004at ISun et alj 
2008; Waelkens et al. 2009). It is also important to empha- 
size that if this emission is local, then the synchrotron halo of 
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FIG. 1. — The 408 MHz synchrotron brightness temperature predicted by 
the wind model presented in Paper I. The predicted synchrotron brightness 
significantly overpredicts the observed emission by approximately a factor of 
two. 

the Milky Way, as inferred from \Beuermann et al.\ ( I79&5I) . is 
substantially overestimated (W. Reich, personal communica- 
tion). 

3.6. Synchrotron Emission from the Wind Model of Paper I 

With all of these components added in, we calculate the 
absolute intensity of synchrotron emission at 408 MHz with 
the original wind model (as presented in Paper I) plus back- 
ground and foreground emission. We find that that model 
substantially over-produces synchrotron emission above the 
plane, as shown in Figure Q] The wind does not dominate 
the synchrotron emission in the disk, but, crucially, falls off 
too slowly with height, violating observational constraints 
strongly near b ~ 6°. As such, the original model in Paper 
I is not in agreement with the observations, and must be dis- 
carded as a "best-fit" to the combined observations. Clearly, 
either the assumed magnetic field is too high, the wind volume 
is too large, or the cosmic -ray electron density is too high to 
fit the observed synchrotron emission. 

We now ask: is there a simple way for the model to be 
brought into closer agreement with both the soft X-ray and 
synchrotron observations, while at the same time, improving 
the basic assumptions of the model of Paper I? 

4. CAN AN IMPROVED WIND MODEL FIT BOTH 
X-RAY AND SYNCHROTRON OBSERVATIONS? 

We now consider two physically-motivated possibilities for 
updating the previously assumed launch parameters of the 
wind as well as one modification for making more robust 
comparisons between the model and X-ray data. These mod- 
ifications are motivated by the overprediction of synchrotron 
emission in the original model but also incorporate current 
thinking about the structure of the inner Galaxy. We then 
check whether such changes may help bring a cosmic-ray and 
thermally-driven wind model in agreement with the observa- 
tions, and consider the effect of these modifications on wind 
driving. 

4. 1 . A Wind From the "Molecular Ring " 



As mentioned in the introduction, Paper I had included the 
assumption that the wind was launched from a range of Galac- 
tocentric radii, from 1.5 to 4.5 kpc. The inner radius was the 
minimum radius at which a wind can be launched in the grav- 
itational potential of the inner Galaxy. The observed edge of 
the soft X-ray emission at / ~ ±30° set the outer radius of the 
wind. In the course of that previous work, it was suggested by 
R. Benjamin (personal communication) that the wind might 
more naturally be consid ered to be launche d from the "5 kpc 
Molecular Ring" (e.g., Jacks on et al.l 12006). which is much 
thinner, with a range of Galactocentric radii of approximately 
AR ~ 1 to 1.5 kpc. This "5 kpc Molecular Ring" is home to 
the largest level of star formation in our Galaxy, and so is a 
reasonable location for a wind launched by hot gas and cos- 
mic rays. So, in this work, we attempt to fit the X-ray emis- 
sion with a wind that covers only 1 kpc on the Galactic disk, 
from R = 3.5 kpc to 4.5 kpc. (This thinner wind may hint at a 
more massive central Galactic bulge; it also increases the im- 
portance of any X-ray emission from static gas in the bulge, 
which we do not include in the model here.) 

For future development of the model, we note that the 
"Molecular Ring" is, of course, not observed to be uniform in 
azimuth, and so an even more realistic assumption for launch- 
ing would include a covering fraction of the wind within this 
annulus. This would require, however, development of a more 
detailed model of the lateral force balance as a function of 
height off the surface of the disk, which is outside the scope 
of this paper. 

4.2. Lowering the Magnetic Field Strength 

In the wind models presented in Paper I, we assumed the 
full magnetic field strength at R = 3.5 kpc of 7.8/iG (F01). 
It is certainly not clear that the entire inferred field strength 
(in the plane of the sky) would be present at the base of the 
wind, and evolve with the wind to high latitudes. In addition, 
if that magnetic field strength were lower, the synchrotron 
emission would of course also decrease, as _B ( t +i V 2 ^ gi.75 
for 7 = 2.5. To test this in the present work, we have de- 
creased the magnetic field strength to 5.2/xG, or 66% of the 
value used in Paper I; this is similar to the local magnetic-field 
strength. We are therefore assuming that only this portion of 
the synchrotron-derived magnetic field at R = 3.5 kpc is verti- 
cal field within the wind, and evolving with the wind to mid- 
and high-latitudes. 

In passing, we point out that the magnetic field strengths 
presented in F01 are, in one sense, lower-limits to the field 
strength. We have found we are able to duplicate the results 
of the emissivity calculation in F01 by using the formulae 
for homogeneous fields perpen dicular to the line of sight (see 
Ginzburg & Svrovatskii| [l965). so that the calculated 7.8fiG 
magnetic field does not include the component of the field 
along the line of sight. On the other hand, it is important to 
note that the emissivity used in F01 is actually the lowest of 
the estimates of local emissivity in Beuerm ann et al.l (11985b 
(as noted in F01), so it is not straightforward to estimate the 
uncertainty in the magnetic field towards the inner Galaxy. 

In addition, it may be possible that the field in the wind is 
not strictly vertical, but rather helical (P. Biermann, private 
communication), as it may reasonably have a strong toroidal 
component due to shear in the Galactic disk. This would 
change the field strength in the plane of the sky, and thereby 
affect the synchrotron emission as a function of Galactic lon- 
gitude. For simplicity, we neglect this effect in the current 
work. 
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4.3. Limiting the Effects of Uncertain Absorption 

In the original fits in Paper I, we found the best-fit to the 
X-ray emission by calculating a x 2 value throughout the en- 
tire latitude range of b = -0° to -90°. (We do not fit the emis- 
sion in the Northern Galactic Hemisphere because of the more 
complex absorption and local emission sources, as discussed 
in Paper I.) However, the fit to the X-ray emission is strongly 
affected by absorption for b > -10° and is dominated by the 
constant background components at b < -40°. Intervening 
absorption was, and is, included in the models, but the dis- 
tance to those absorbers is unknown, and so the absorption 
is calculated as a screen of material located between the ob- 
server and wind. This is perhaps a reasonable first assump- 
tion, but we can minimize the impact of that uncertainty in the 
fit by only calculating the x 2 value between -40° < b < -10°, 
which we do for all of the following models. This also limits 
the impact of uncertainty in the stellar contribution to X-ray 
emiss ion at low-latitudes (Masui et al. 2009; Revnivtsev et al.1 
2009), although such emission is more significant at higher 
energies. (In practice, we find that these modifications do not 
significantly change the parameters of the best-fit models in 
the previous or current work.) 

4.4. Putting it All Together: an Improved Wind Model 

Decreasing the wind model's footprint, so that it stretches 
only from 3.5 to 4.5kpc on the disk, and attempting a smaller 
magnetic-field strength (B = 5.2pG) for the wind, we inte- 
grate the hydrodynamic equations for a range of initial (at 
the base of the wind) thermal-gas pressures, thermal gas den- 
sities, and Zbreak values to find the best-fit model (using the 
same methods of Paper I; see the Appendix) to reproduce the 
X-ray emission. This search entailed two large-scale searches 
through parameters space that integrated the hydrodynamic 
equations of motion for 13,475 models each. The first search 
spanned two orders of magnitude in thermal-gas density at the 
base of the wind, over one order of magnitude in thermal-gas 
pressures at the base of the wind, and Zbreak values spanning 
a factor of two (in 35, 35, and 1 1 steps, respectively), cen- 
tered roughly near the best fit parameters from the wind model 
in Paper I. The best-fit to the soft X-ray emission from this 
scan was then input to a parameter-search of another 13,475 
models distributed over a slightly smaller range in parameters 
(thermal-gas density ranging over a factor of ten, thermal-gas 
pressure ranging over a factor of 2.25, and Zbreak ranging over 
a factor about 1.5), to zoom in near the best fit. A final level 
of refinement in the parameter search was run to scan 1331 
models in those three parameters (only varying each parame- 
ter by 40%, 20% and 20% around the previous best-fit value 
for n gi o, Pgfl, and Zbreak, respectively). 

The new best-fit model, found by this procedure, is shown 
in Figures |2] and [3] with the parameters for the fit given in 
TableQ] We also show the velocity and density profiles of this 
wind in Figures [4] and 

This model still improves on the static model of Alm v et all 
(2000) with respect to \ 2 - First, calculating x 2 over the en_ 
tire range of longitude (-90° < b < 0°, which is not the re- 
stricted region based on which the model was selected), the 
new model improves on the static polytrope's \ 2 by a factor 
of 1.6 in the ROSAT 'R4' band and 2.1 in the 'R5' band. It 
is important to note that this fit is not as good as the origi- 
nal fit from Paper I, however: the new model overpredicts the 
emission in the X-ray 'R4' band by 15% at intermediate lat- 
itudes (b ~ -30°), and slightly overpredicts the emission in 



'R5' at low latitude (I — -10°fl (We note that this best-fit 
model is defined as the model with lowest x 2 , calcu lated in 
the region from -40° < b < -10° as defined in j4~3l but this 
only changes the best-fit model parameters by < 3% from the 
fit calculated with x 2 defined over the entire latitude range.) 

How have the changes in AR and B affected the wind 
model? In order to fit the X-ray emission while being 
launched from a thinner annulus, the gas density must in- 
crease at the base of the wind (here, by ~ 30%). The gas 
pressure must increase as well, to keep the wind at approx- 
imately the same temperature (which is constrained by the 
relative emission in the 'R4' and 'R5' bands). The gas pres- 
sure has to increase still further for another reason, as well: 
as the magnetic-field strength has decreased and the density 
has increased, the power that is channeled from the cosmic 
rays to the gas also decreases (the power goes as Va • VP C1 oc 
■ VP C1 ). In the model of Paper I, energy from cosmic-rays 

played a significant role in heating the gas at mid-latitudes 
of b ru 15°. But, as B decreases and p increases, the power 
transmitted from the cosmic -ray component decreases, and 
the base gas tem pera ture must increase still further to com- 
pensate (see Eq. IA3l of the Appendix). This leads to an in- 
crease in the gas pressure at the base of the wind of approxi- 
mately 104%. These are the dominant changes to the model 
as a result of the synchrotron constraints on the wind model. 

Turning to other parameters of the model, the lengthscale 
for flow tube expansion (zbreak) has decreased by ~ 25% (this 
parameter sets the scale-height for the expansion of the as- 
sumed flow tubes, and so, in a sense, represents an impor- 
tant set of assumptions about the magnetic geometry of the 
Galactic halo, which helps to define the acceleration profile). 
Meanwhile, the cosmic-ray pressure at the Galactic midplane 
has remained constant, by assumption. To further compare 
models, the mass outflow rate in this wind is 2.2M /yr (very 
similar to the 2.1M /yr for the previous model), and the ter- 
minal velocity is ~ 570km s -1 , significantly lower than the 
760 km s" 1 for the wind of Paper I. Overall, the model is some- 
what similar to the past one in its ability to fit the X-ray emis- 
sion, and in its parameter values, with the notable exception 
of strongly increased gas pressure. It is also important to point 
out that the initial velocity of gas in this simple wind model 
is 252 km s" 1 ; such a high initial velocity is the result of our 
assumption that all of the energy is input as a delta function at 
the Galactic midplane, and so strongly indicates the need for a 
more detailed model of more gradual momentum and energy 
input to the wind. 

Given this best-fit model to the soft X-ray emission, we can 
compute the synchrotron emission for the model. We com- 
pare the resultant brightness temperature to 408-MHz obser- 
vations in Figure|6] This figure shows the averaged brightness 
temperature from the 408 -MHz survey compared to the wind 
model outlined above (when all background and foreground 
sources have been added); both the data points and model 
curve are the result of an average over -30° < I < 30°. The 
wind model now successfully fits the observed synchrotron 
emission at the mid-latitudes of -25° < b < -10°. 

5 It is important to stress that this assumes that the wind model plus back- 
grounds explains all of the X-ray emission observed towards the center of 
our Galaxy; it is very conceivable (especially now that the wind is much thin- 
ner than in Paper I) that other components may help explain the emission 
(such as the loca l ISM, the Galactic bulge, or even a bulge-related wind, see 
ITang et aU2009l) , and therefore the requirements for the wind could be some- 
what relaxed, and the strength of emission in the wind therefore decreased. 
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FIG. 2. — The best- fit model for a Galactic wind that stretches only between 
3.5 and 4.5 kpc in Galactocentric radius with B = 5.2/xG; this comparison is 
for the ROSAT 'R4' band, centered roughly on 0.65 keV. The x 2 for this fit is 
2335 for 87 degrees of freedom, yielding a xt of 26.8; this is slightly worse 
but comparable to the value of x 2 , = 19.0 from the best-fit R4 model in Paper 
I. 



We point out that this model fits simply due to the decreased 
AR and B; as such, the model presented here is not unique: 
we have essentially modified the model to include more phys- 
ically motivated launching conditions, finding a model that 
satisfies the observational requirements. There are, of course, 
a range of models that might also fit the observations, given 
corresponding changes in B and P c the range of possible 
models in this parameter space are shown in Figure [7] This 
figure displays the range of parameter values for P c o and B 
for which winds are successfully launched; the color contours 
show the region of escaping winds with the colors encoding 
the x 2 value for the fit to the two bands of soft X-ray emis- 
sion, while the red lines show an estimate of the factor by 
which models in that parameter space would over- or under- 
predict the synchrotron emission: the red lines show the in- 
crease in synchrotron emission by a simple scaling of the best- 
fit synchrotron emission by n cr e • B 115 . This figure shows the 
range of magnetic-field strengths and cosmic-ray pressures 
that would approximately satisfy the synchrotron (with the 
red line labeled '1.0') and X-ray observations. The figure 
is also interesting in its own right, as it shows the wide range 
of magnetic-field strengths that can help to launch a wind. 
The biggest difference between the various models over the 
range of B is that the models with larger B can more easily 
reproduce the large scale-height of soft X-ray emission due to 
distributed heating by the cosmic-ray protons. 

Note that we have added a constant foreground component 
(as well a s oth er already known background sources, men- 
tioned in jj3. 41 to our model in Figure [6] If the high-latitude 
emission is not local, the constraints on AR and B would be 
relaxed, and larger AR and B would be allowed. 

This fit assumes that all other parameters are constant, to 
attempt to limit the number of variables. In particular, the 
cosmic -ray pressure at the base of the wind, P c .o, has not been 
allowed to vary in any fits (although Fig.|7]displays the param- 
eter space of P c o and Bo to show how joint variations affect 
synchrotron observations). In Figure [8] we show the param- 
eter space where P c ,o and P gi o are allowed to vary together; 



FIG. 3. — As in Figure [2] but for the ROSAT 'R5' band that is centered 
roughly on 0.85 keV. The \ 2 f° r this fit is 4522 with, again, 87 degrees of 
freedom; xi i- s then 52.0, which again is only slightly worse fit than the 
xl = 48.9 fit from Paper I. 
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FIG. 4. — Velocity vs height for the best-fit model (parameters given in 
Tableff}. The solid line shows the thermal-gas velocity within the wind, the 
dashed, red line shows the composite sound speed as a function of height, and 
the dot-dashed, green line shows the variation in Alfven speed with height. 
The increase of velocity (even through the critical point, where v = c*) is 
fairly standard, although this wind exhibits a slight decrease in velocity at 
large distances, as the increased wind radius encompasses more and more of 
the Galaxy's mass. 



TABLE 1 
Wind Parameters: New and Old 





New 


Old 




Parameter 


Value 


Value 


Fixed? 


Pg.o/ks 


4.0 X 10 4 K cm" 3 


2.0 X 10 4 K cm" 3 


Varied 


Po 


9.0 X 10~ 27 gem -3 


7.1 X l(r 27 gem" 3 


Varied 


£break 


4.0 kpc 


5. 2 kpc 


Varied 


Rq Range [Galactocentric] 


3.5 to 4.5kpc 


1.5 to 4.5 kpc 


Fixed 


Pcfi/kB 


2.2 X 10 4 KcnT 3 


2.2 X 10 4 KcnT 3 


Fixed" 


Bo 


5.2 ^G 


7.8 ^G 


Fixed" 


a 


2.0 


2.0 


Fixed 



11 Varied in FigurefTJonly, and not fit to the data. 
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FIG. 5. — The temperature and density in the thermal gas as a function 
of height in the wind. The temperature is given by the solid, black line; as 
expected, as the wind accelerates and expands into the assumed fiowtube, 
the temperature decreases. The temperature does not decrease as quickly as 
without cosmic-ray heating, however. Also in this plot, we show the particle- 
number density as a function of height with the dashed, green line. 



Averaged Temp over lil < 30 Degrees, S. Hemi. 



Averaged Model Synchrotron Emission + Bkgds 

O Hoslem et ol. (1982) 408 MHz Doto 

Model Synchrotron Emission 

CMB Background 

Extragalactic Background 

Isotropic Background 




Galactic Latitude [degrees] 



FIG. 6. — The integrated 408 MHz synchrotron brightness temperature 
given by a wind that covers a smaller area on the Galactic disk than the model 
in Paper I; this wind extends only from 3.5 kpc to 4.5 kpc on the disk, and has 
B = 5.2/xG. 



this plot shows the importance of adding cosmic -ray pressure: 
the area of allowed winds (shown by the color contours) in- 
creases markedly for P c fl > 10 4 cm" 3 K. In that regime, Galac- 
tic winds exist for lower thermal pressures than would other- 
wise be allowed. This plot also shows the best fit (from Ta- 
ble [TJ; however, interestingly, when P c o is allowed to 'float', 
a fit is preferred where P c o and P g o are approximately equal. 
However, this results in a cosmic -ray electron density that, 
given our stated assumptions abou t the ratio of the cosmic -ray 
electron-to-proton fraction Q3.lt and other parameters in the 
fit, would violate the radio-synchrotron constrains. Of course, 
if the proton-to-electron ratio was 100:1 instead of 50:1, as 
assumed here, that would relax this constraint on the wind 
greatly. 

Finally, in Figure [9] we show the allowed winds as a func- 



FlG. 7. — The parameter space of the best-fitting model as a function of 
the magnetic-field strength at the base of the wind, Bq, and the cosmic-ray 
pressure at the base of the wind, P c ,o. The color-contour-delineated region 
denotes those winds which can be launched from the Galactic potential, and 
the color contours show the combined, total \ 2 f° r fitting the X-ray emission 
in the ROSAT 'R4' and 'R5' bands. The labeled red lines indicate the trend 
of synchrotron emission, with the number giving the estimated level of syn- 
chrotron emission above that observed in the best-fitting model. One can see 
a range of models with radio emission of the level of that observed (the red 
line marked '1.0') with relatively low x 2 values. The best-fit model to the soft 
X-ray observations (with P c u fixed) is indicated by the green 'x'. NB: The 
red lines here give the synchrotron scaling of /„ oc « C r,eB'' 75 appropriate for 
7 = 2.5, and do not represent the results of a set of synchrotron simulations, 
which would yield modified curves due to variations in wind acceleration 
profiles impacting the run of density with height. For winds in this survey, 
the parameters P„o, no, Zbrcak and a were held fixed. 

tion of the wind parameters no and P g Q. This figure is much 
like the corresponding Figure 10 in Paper I, albeit with a 
'slimmer' region of allowed winds due to the increased im- 
portance of gas pressure vs. cosmic -ray pressure in the newer 
model (due to the smaller magnetic-field strength). 

4.5. Another Diagnostic: the Wind's Brightness 
Temperature Spectral Index 

There is another possible constraint on the wind from syn- 
chrotron observations: the temperature spectral index. Af- 
ter much remarkable work on surv ey and survey calibrations 
over the last few decades (e.g., Hasla m et alj 119821 : iReich] 
[1981 iReich & Reich][T986t IReich et alj|2001l) maps exist of 
the syn c hrotron spectral in dex over the sky (IReich & Reich! 
Il988alfbl : IReich et al.|[2004b . In preparation for a new analy- 
sis of the spectral index (s ome early results have already ap- 
peared in Reich et al. 2003), we can make predictions of the 
temperature spectral indices that would be produced by the 
wind. We calculate the brightness temperature spectral index 
from predicted 408 MHz and 1420MHz. This spectral index 
is defined as: 

where for this analysis, v x = 408 MHz and v 2 = 1420 MHz. 
As in the case of total intensity, we calculate this tempera- 
ture spectral index for the average emission in the model and 
in the data for 1 = 0°, and over the range of latitude where 
the wind emission seems to dominate the observed emis- 
sion, 10° < \b\ < 25° (and also, where d iscrete sources are 
much less numerous). For a similar region. lReich et alj d2004l) 
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FIG. 8. — The parameter space of successful outflows as a function of gas pressure at the base of the wind, P„ q, and cosmic-ray pressure at the base of the wind, 
P c q. The color contours outline the region of allowed winds, with the color code specifying the mass outflow rate in the wind. However, in this figure, the black 
lines denote the total pressure in the wind (cosmic-ray pressure plus thermal pressure). As expected, the mass outflow rate generally increases with total pressure 
in the wind. This figure also shows the important role played by cosmic-ray pressure, though: cosmic-ray pressure above P z o ~ 10 5 cnT 3 K result in successful 
winds for values of P s q that would ordinarily not allow outflows. For winds in this survey, the parameters tig, Bq, Zbreak an d a were held fixed. 



found the temperature spectral index /3 ~ 2.65 at b ~ 30° (see 
their Fig. 3; this measurement includes all components along 
the line of sight). 

The temperature spectral index for the best-fit wind model 
(alone) is shown as the solid line in Figure[l0] When we con- 
sider the spectral index of the wind (without the foreground 
emission component accounted for), the wind spectral index 
gradually steepens, as expected for the cooling cosmic-ray 
electron population: the wind model has (3 ~ 2.85 at b = 10°, 
rising to (3 ~ 3 . 1 by b = 40 °. This is significan tly steeper than 
the spectral indices seen in lReich et al. I (120041). which approx- 
imately span the range of 3.0 to 2.6, falling with height above 
the disk. 

This source of "flattening" of the observed synchrotron 
spectral index has been the source of some d ebate (e.g., 
iReich & Reichhl988at|Pohl & Schlickeiserll 19901) . In this de- 
bate, we would like to emphasize the importance of multi- 
ple emission compon ents that may lie along the line of sight 
(as also mentioned in Pohl & Schlickeiserl l 1 9901) . Our model 
(when considered alone) clearly shows a steepening of the 
spectral index with height, but when we calculate the spec- 
tral index of the wind plus a constant (again, perhaps local) 
foreground, we find a nearly constant spectral index; this is 
shown by the dashed line in Figure [10] We note that this is 
the simplest foreground possible, but it is sufficient to signifi- 
cantly change the run of spectral index with height. 

As such, until we have a more detailed understanding of 
background and foreground synchrotron sources, it will be 
difficult to draw definitive conclusions about any one com- 



ponent. What is of course needed here is an improved un- 
derstanding of the nature of the various components and how 
they work together, in order to constrain the various physical 
processes and dynamics of our Galaxy ( Strong et al. 2004a; 

ICoxll2005T) . 

5. RESULTS & DISCUSSION 

We have continued to build on the results of previous 
work, and on our own Paper I, by applying radio-synchrotron 
constraints derived from the 408 MHz all-sky survey of 
Haslam et al. ( 1982) to our cosmic -ray and thermal-gas pres- 
sure driven wind model. The original model of Paper I 
overpredicts mid-latitude synchrotron radiation (although it 
would approximately reproduce the observations if a cosmic- 
ray proton-to-electron ratio of 100:1 was adopted; we retain 
the locally observed ratio of approximately 50: 1 to place the 
most stringent constraint on the wind model). However, a 
wind model can be found that is consistent with the observa- 
tions and with more physically-plausible launch conditions, 
such as a thinner wind with AR ~ 1 kpc and with smaller 
B ~ 5^G. 

Both the trend of star formation with Galactocentric ra- 
dius and radio-synchrotron observations point to a wind that 
is launched from a smaller "footprint" than hypothesized in 
Paper I. This range of radii seems plausible if we envision the 
wind as launched from within the inner spiral arm, or the inner 
molecular ring, which has been observed in this radius range. 
We also found that a smaller magnetic field strength was re- 
quired: the vertical field in the wind can be decreased from 
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FIG. 9. — The parameter space of allowed models as a function of gas pressure at the base of the wind, P„ o and density, hq. The color contours show the range 
of models where plasma escapes the Galaxy, with the color corresponding to the mass outflow rate, given in the color bar on the right. The solid lines are marked 
with the gas temperature along that line; the red-dashed line give the lower (7j ow ) and upper (Thigh) limits for the gas temperature of a purely thermally-driven 
Galactic wind. The best fit is shown with the yellow 'x', surrounded by the region where \ 2 < 2 ■ X^j n > giving an indication of the region of minimally acceptable 
fits to the X-ray emission. Compared to the similar figure in Paper I (that paper's Fig. 10), the region of allowed winds is smaller at higher pressure due to the 
relative dominance of gas pressure in this winds. For winds in this survey, the parameters P c .o, Bq, zt,reak and a were held fixed. 
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FIG. 10. — The temperature spectral index, /3, for the Galactic wind model 
along / = 0° . The solid line shows the trace of the spectral index for the wind 
model alone, while the dashed line gives the spectral index for the wind with 
a foreground emission component. 



7.8/iG to 5.2/iG and fit the observations. The resulting wind 
model predicts a level of synchrotron emission that agrees 
with the level observed at high latitude, does not overpredict 
synchrotron emission at low- to mid-latitudes (b ~ -15°), and 
reproduces the X-ray emission in the ROSAT 'R4' and 'R5' 
bands. However, the fit in the 'R4' band, while still improv- 



ing on previous static models, does not fit as well as model in 
Paper I. We have suggested that this may be due to the lack, in 
this model, of X-ray emission from static gas. Still, the con- 
straint of the synchrotron emission requires that the wind have 
a smaller radial extent, and have a smaller magnetic field; the 
'best-fit' model, in this paper, represents the best agreement 
between these two different constraints (while keeping P c q 
constant). 

While comparing the synchrotron observations to the wind 
model, we have found that a significant fraction of high- 
latitude (b < -40°) synchrotron radiation is consistent with 
a local source rather than Galactic - halo emission. T his is in 
agreement with other work (see, e.g. lSun et al.ll2008l and ref- 
erences therein), and hin ts that the synchrotron halo model of 
iBeuerm ann et all d 1985b is perhaps a signific ant over-estimate 
of the actual Galactic halo. As the model of IBeuermann et all 
d 1985b is still in widespread use, this seems an important point 
to stress. We include this local emission to again strongly 
constrain this Galactic wind model, to see if it can survive. It 
appears that it can. 

Choosing this wind model has other implications for the 
effect of a wind on the Galactic central region. We wish to 
point out, however, that we do not consider the diffusion of 
electrons throughout the Galaxy; only the vertical advection 
of cosmic rays within a small range of Galactocentric radii 
is considered. Given this small range of launching radii, and 
since the radial diffusion rate of cosmic-rays is much less than 
the vertical advection rate, the present wind model would not 
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advect cosmic rays from a widespread range of Galactocen- 
tric radii. Indeed, we envision that, since cosmic rays stream 
mostly along magnetic field lines, this wind's advection of 
cosmic rays would not impact the cosmic -ray distribution in 
the solar neighborhood. 

This advection of cosmic rays from the Galactic mid- 
plane could have other effects, though: in Paper I, it was 
hypothesized that such a wind could help evacuate cosmic 
rays from the central Galactic disk where 7-ray observations 
seem to indicate that cosmic -ray protons are not as abun- 
dant as the apparent su pernova rate i n that region of the 
Galaxy would sugg est (Bloemenl [l989t iBloemen et al.lll993t 
iBreitschwerdt et all 120021) . This type of wind and its effect 
on the cosmic-ray population ha s been studied in more de- 
tail by Gebauer & de Boerl d2009l) . The relative dearth of cos- 
mic rays towards the center of the Galaxy might then require 
an additional component, such as a sig nificant change in th e 
W C o-to-N(H 2 ) factor, as suggested by Strong et al. (2004b). 
This model will also require more localized supernova power 
to launch the wind. In the wind model presented in Paper 
I, the wind required approximately 2.1 times the supernova 
power that is observationally inferred in the radial range be- 
neath it (R = 1.5 to 4.5 kpc), using the estimate of the super- 
nova rate as a function of R in F01. Using the same method, 
the wind model presented in this paper requires approximately 
2.7 times the supernova power in the AR = 1 kpc range of radii 
that it is launched from, assuming that each supernova pro- 
duces 10 51 ergs~ 1 (using instead the supernova rates vs R of 
iMcKee & Williamsll 19971 the wind requires 1 .9 times the es- 
timated supernova power). While, as mentioned in Paper I, 
thermal conduction and clumping of the outflowing gas may 
help reduce this estimate, this calculation shows that the syn- 
chrotron data has perhaps pushed the simple wind model to 
the limit of its applicability to the Galaxy. However, this es- 
timate is also very dependent on the radial distribution of su- 
pernovae; if uncertainties in the distances to supernovae are 
artificially spreading out the distribution of supernovae, then 



the power in the Molecular Ring may be underestimated. All 
of these considerations point to the need for not only a more 
complete model but also for continued comparisons of the 
model with other observations. 

Overall, this paper has shown that the wind model can sat- 
isfy both the synchrotron observations and can be made more 
compatible with what is already known of Galactic structure. 
This work satisfies those observations and constraints by lim- 
iting the radial extent of the wind and magnetic field strength 
within the wind. The most direct way to continue testing the 
model is to produce 7-ray predictions, and compare those pre- 
dictions with Fermi/LAil observations of the diffuse Galactic 
7-ray emission. 
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APPENDIX 



WIND MODEL EQUATIONS 

For completeness, we review here the equations for the cosmic-ray and thermal-gas pressure driven wind from lEverett et al.l 
(2008); please refer to that paper for further details. 

Our wind model builds on the equations first presented by IBreitschwerdt et al.l ( 1199 11) and Breitschwer dt et al.l (TT993). This 
model is a ID, semi-analytic model that treats both the thermal gas and cosmic -ray components as fluids. The development of 
the outflow with height is governed by the equation of mass conservation (assuming no mass is added as the wind flows out of 
the Galactic plane): 

^-(pvA) = (Al) 
as 

where p is the gas mass density, z is the height above the Galactic plane, and A is the cross sectional area of the wind. Since this 
is a ID system of equations, we must prescribe the cross-sectional area as a function of height: 



A(z)=A Q 



1 + 



£break 



(A2) 



where Ao is cross-sectional area of the flowtube at the Galactic midplane, and a is the power-law governing the divergence of the 
flowtube; the flowtube has a roughly constant cross section A ~ Ao until approximately the height Zbreak, where the area increases 
as z a ■ We chose a = 2 to mimic a spherical divergence with height above Zbreak and leave Zbreak as a parameter in our fitting 
process, as its value is not known beforehand. 
The thermal-gas and cosmic -ray pressure change with height via the relations: 

\ 

(A3) 
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Everett, Schiller & Zweibel 



dP c lc P c M A +\dp 

and = — , (A4) 

dz p M A +l dz 

where P g and P c are the thermal-gas and cosmic-ray pressures, c g is the speed of sound, 7 g and -f c are the polytropic indices for 
the thermal gas and cosmic rays (set to 5/3 and 4/3, respectively), and Ma = v/va is the Alfven Mach number for the gas (v is 
the velocity of the gas and va is the Alfven speed). 

The magnetic-field strength evolves with height by requiring that the magnetic flux remains constant; hence, B(z)A(z) = 
[constant]. We note that the magnetic field is important here not as a direct source of pressure or tension, but as a "conduit" 
of sorts, through which the cosmic rays communicate momentum and energy to the thermal gas. This is why we require a vertical 
magnetic field; a magnetic -field component must lie along the direction of the cosmic-ray pr essure gradient in order to excite the 
streaming instability (e.g., Wentzel 1968; Kulsrud & Pearcdll969fclBreitschwerdt et al.lfl99 lb . 

The terms in these equations (Eq. |A3| and lA4l ) can be understood as follows. The first term in Equation IA3I de scrib es how 
the gas pressure changes due to expansion and acceleration of the gas in the fiowtube. The second term in Equation |A3| couples 
the cosmic -ray pressure to the gas pressure; this term represents the damping of cosmic-ray generated Alfven waves which then 
heat the gas and help drive the wind. This type of immediate wave damping is important as rapid wave-damping mechanisms 
are known to occur in the ISM, and this also limits the growth of the Alfven waves such that the perturbation in magnetic-field 
strength remains below the background field strengt h. This idea was firs t intro duced (but not yet wid ely used, since the Alfven 
wa ve growth was small below the critical point) in Breitschw erdt et alJ (1199 ll) and was al so u sed in Zirakas hvili et alJ (119961) 
and Ptuski n et al.1 (119971) ; we discuss this in more detail in §2 of Paper I. Finally, Equation lA4l describes the loss in cosmic-ray 
pressure due to adiabatic expansion, wave momentum transport to the gas, and to wave generation by the cosmic rays. 

The velocity of gas in the wind is the derived by solving the wind equation, which we write as: 

dv 2 dp 

p v -r +c *-r=-pg- ( A5 ) 

dz dz 

This equation of motion is much like the equation of motion for models that rely only on thermal-gas pressu re, except that instead 
of a gas sound speed, c g , we use c*, which gives a "composite sound speed" (see lBreitschwerdt et al.llT99 1 ). defined as: 

, d(P s + P c ) 

cl = s , (A6) 
dp 

where P g and P c were given in Equations I A3 1 and I A4I 

To solve for the functions v(z), p(z), Pg(z), P c (z) and B(z), we must integrate the above equations from the Galactic midplane 
to a distance far above the galactic midplane, correctly threading the critical point that occurs when v = c». When starting 
the integrations from the midplane, we set (from other models and observations) and fix the value of the midplane cosmic-ray 
pressure, the magnetic field strength, and the power-law a that governs the opening of the fiowtube cross-sectional area (see 
TableQ]for values). We do not have strong constraints on the values of the midplane thermal-gas pressure, the hot-gas midplane 
density, and the height where the fiowtube starts to rapidly expand (ztreak) so we leave those as free parameters which are then 
constrained by the soft X-ray and radio observations. These parameters are also considered in some detail in Paper I and in §0 
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